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ABSTRACT 

We present a system of X-ray photometry for the Chandra satellite. X-ray photometry can be a 
powerful tool to obtain flux estimates, hardness ratios, and colors unbiased by assumptions about 
spectral shape and independent of temporal and spatial changes in instrument characteristics. The 
system we have developed relies on our knowledge of effective area and the energy-to-channel conver- 
sion to construct filters similar to photometric filters in the optical bandpass. We show that the filters 
are well behaved functions of energy and that this X-ray photometric system is able to reconstruct 
fluxes to within about 20%, without color corrections, for non-pathological spectra. Even in the worst 
cases it is better than 50%. Our method also treats errors in a consistent manner, both statistical as 
well as systematic. 

Subject headings: X-rays: general — techniques: photometric 



1. INTRODUCTION 

Photometry is one of the most widely used, relatively 
simple, tools used in describing and categorizing astro- 
nomic al objects. Standardization by Jo hnson fc Morganl 
(1953) and subseq uent ad di tions in the optical and in- 
frared (see e.g. iBesselll ()2005l )) have allowed com- 
parisons between measurements by different telescopes 
and instruments without bias. Important applica- 
tions of o ptical photometry in cl ude s tellar classifications 
(see e.g. Uohnson fc Morganl (|1953l )l. galaxy redshifts 
(jPuschell et al.lll981l ), an d the discovery o f the most dis- 
tant quasars in the SDSS (|Fan et al.l fl999) . In particular 
photometry is important for sources too faint to extract 
detailed spectra, i.e. most sources, given the almost uni- 
versal increase of source numbers to faint fluxes. 

Although optical astronomy is a much older discipline 
than X-ray astronomy, optical photometry was estab- 
lished only about 20 years before the beginning of X-ray 
astronomy in the f 960s. Given the enormous success of 
optical photometry it seems obvious to use it as a prece- 
dent and try to duplicate its success in other wavebands 
beyond UV and infrared. The X-ray band can be defined 
as reaching from about 0.1 keV to a few hundred keV, 
spanning almost 4 decades of frequency - although most 
work concentrates on the 0.2-20 keV band - compared 
with 2 octaves in the optical. Moreover, an important 
difference between optical and X-ray is the typically low 
number of photons in X-ray astronomy. The X-ray range 
is photon starved such that sources with a few hundred 
counts are considered bright in X-rays. This limitation 
increases the importance of broad-band photometry in 
the X-ray band. 

While X-ray astronomy has used relative, mission- 
specific photometry for most of its existence, there is, as 
yet, no standard X-ray photometric system. The use- 
fulness of an X-ray photometric system is evident al- 
ready from the use of these somewhat idiosyncratic en- 
ergy bands. The bands used in the past have been chosen 
for specific purposes, e.g. to use color -color diagrams to 
diagnose X-ray binary spectral states (|White fe; Marshall 
11984 lHasinger fc van der Klisl I1989D where, in the lat- 
ter, the energy bands are different for each source. Even 



so, the resulting color-color diagrams have immensely in- 
creased our knowledge of X-ray binary spectral /accretion 
state s (e.g. iPrestwich et al.l 120031 : iGierlihski fc Newton! 
2006). Thus a standard X-ray photometric system is 
highly desirable in X-ray astronomy in order to cross- 
compare observations of the hun dreds of thous ands of 
sour ces being cataloged by XMM (|Watsonl 120071 ) . Chan- 
dra (|Fabbiano et alj |2007| ). and other missions. Even 
within a given mission different types of CCDs (XMM) 
or changes in operating temperature, gain or contamina- 
tion {Chandra ) mean that simple count rates cannot be 
used. 

However, there have been complicating factors in es- 
tablishing photometric energy bands beyond individual 
observations. One cause of this lack of standardiza- 
tion is that the energy ranges covered by different X-ray 
satellites and instruments differ widely. For example, 
RXTE/PCA, G1NGA and EXOSAT bands have practi- 
cally no overlap with ROSAT bands; and ASCA, XMM 
and Chandra bands are somewhere in the middle. Fig. [T] 
shows a selection of energy bands used by different au- 
thors for different X-ray satellites 1 . Most X-ray missions 
with focusing optics cover the energy range from ~0.1 
keV to 10 keV. 

Another reason for the lack of a standard photomet- 
ric system is that in X-rays there are no bright constant 
point sources in the way that stars can be used for cali- 
bration like in the optical. 

But the most fundamental cause for the lack of an X- 
ray photometric system has been the limited spectral re- 
solving power (K=E/AE ~1) of proportional counters, 
which were used in X-ray astronomy from the earliest 
days through to ROSAT and RXTE. A resolution of R~l 
allows no clean separation of energy bands, and different 
spectra with similar flux will give widely different flux es- 
timates in any chosen band. In optical terms, the "color 
correction" is very large. However, with the introduction 
of X-ray CCDs in ASCA (jBurke et al.lfl99"l this limita- 

1 For more information on X-ray satellites, see the HEASARC 
web page http://heasarc.gsfc.nasa.gov/docs/observatories.html A 
list of energy bands and references is given in the on-line version 
of this paper. 
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tion has largely gone away. X-ray CCDs have R>10, so 
comparable to the R^6 of broad band optical photome- 
try. It seems thus timely to consider the introduction of 
an X-ray photometric system. Therefore we have inves- 
tigated how good a photometric system can be created 
for Chandra ACIS observations and, by extension, for all 
other X-ray CCDs. We report the encouraging results in 
this paper. 
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Fig. 1. — Energy bands used in various publications for X-ray 
satellites. The different energy bands for different satellites are due 
to energy coverage of the instruments. The colors represent the 
width of the soft/medium/hard bands used in the corresponding 
papers. A table of energy bands and references are given in the 
on-line version of the paper. 



2. MEASUREMENTS WITH Chandra ACIS 

The ACIS (Advanced CCD Imaging Spectrometer) in- 
strument is one of the two detectors on Chandra . ACIS 
is arranged as an array of 10 1024 by 1024 pixel CCDs, 
one as a two by two array (ACIS-I), and one as a 1 by 
six array (ACIS-S). The instrument allows simultaneous 
high resolution imaging and moderate resolution spec- 
troscopy. Two CCDs on ACIS-S are back-illuminated 
(ACIS-S 1 and ACIS-S3), the rest is front-illuminated. 

2.1. Energy bands 

The black vertical lines in the Fig. [T] are the energy 
bands we have chosen. They are roughly an octave wide, 
though the hard band is somewhat wider. 

TABLE 1 

Standard photometric bands for Chandra 



s 

soft 


Ml 

medium-soft 


M2 

medium-hard 


H 
hard 


0.3-0.5 keV 


0.5-1.0 keV 


1.0-2.1 keV 


2.1-7.5 keV 



Most publications use similar bands as can be seen in 
Fig.ffl 

The low energy bound is chosen to be 0.3 keV because 
the steepness of the ACIS effective area curve for E < 0.3 
keV leads to large dependencies on spectral shape and, 
secondly, the effective area and photon energy-to-channel 
conversion for the S-3 chip are not well calibrated below 
0.23 keV (P. Plucinksy, priv. comm.). The iridium edge 
at ^2.1 keV defines the separation between the second 
medium and hard band. Note that other X-ray mirrors 
are coated with other high Z elements (e.g. Au, Pt) but 
these have edges at similar energies (2.3 keV and 2.13 
keV, respectively) . The high energy boundary is defined 
by the rapidly rising ACIS background and falling effec- 
tive area at energies above ~7.5 keV (Chapter 6.16 in 
ICXCl ([2001 ). 

The exact location of the boundaries will not strongly 
affect the analysis of source colors, except for pathologi- 
cal spectral shapes or strong lines at or close to bound- 
aries. For Chandra a much stronger effect, at least in the 
soft(er) band(s), is the build up of contaminating mate- 
rial on the front of the ACIS, which reduces the effective 
area significantly compared to the case o f no contamina- 
tion ( see O and iPlucinskv etai] (|2003f ) ; iMarshall et~aTl 
((III). 

The number and boundaries of energy bands are of 
course not immutable, optical photometry also has a va- 
riety of more or less different bands for different scientific 
purposes. For example narrow band filters at neutral or 
hydrogen-like iron lines or silicon lines can be of interest. 
Our program to compute correction factors for the count 
to flux conversion is not limited to the above mentioned 
energy bands. However, the use of a common standard 
system across missions is highly desirable as explained in 
SeclU 

2.2. ACIS ARF 

The ARF (Ancillary Response File) describes the ef- 
fective area of a telescope at a given energy. All effects 
related to the probability of detecting a photon with 
the telescope and detector (quantum efficiency, block- 
age, mirror area, vignetting) are combined in the ARF. 
In an ideal system the effective area would be indepen- 
dent of energy. However, the shape of the effective area 
curve of X-ray telescopes varies far more than for opti- 
cal telescopes, with variations of factors of a few, even 
across the octave-wide energy bands used here. Optical 
telescopes also cover a much smaller logarithmic range of 
photon energies. 

For the back-illuminated (BI) ACIS chips of Chandra , 
there is significant time variation in the effective area 
at low energies. Due to resublimation of material on the 
chips, sensitivity at low energies, up to 0.6 keV, decreased 
by a factor of ~30 since the beginning of the mission. The 
sensitivity of the front-illuminated (FI) chips also suffers 
from accumulation of material but the reduction of low 
energy effective area is not as pronoun ced as the area 
is alr eady quite small for these chips (jPlucinskv et al.l 
2003). This is shown in Fig. EH 

For more details about changes in the Ch andra ARF 
see lSchwartz et al.l (|2000t ); IZhao et al.l (|2004| ). 



Because of the similarity of the XMM and Chandra en- 
ergy ranges wes. will use these bands as a starting point. 
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Fig. 2. — Effective are curves for the front-illuminated ACIS-I3 
chip (top) and the back-illuminated ACIS-S3 chip (bottom) at the 
aimpoints for the years 1999 to 2008. Years are shown in different 
lines styles: dotted-1999, dashed-2000, to solid 2008. Years 2004- 
2007 are not shown because they are almost identical to the 2008 
curve. The vertical dotted lines delimit the energy bands. 

The RMF (Redistribution Matrix File) is a matrix 
that redistributes incident photon energies E (raw en- 
ergy channels) to PHA or PI channels or observed energy 
E' . Neglecting the effective area the observed spectrum 
in a channel Sp(E') is related to the true spectrum at a 
given energy E by the convolution of incident spectrum 
with RMF: 



Sp(E') = Sp{E) ® RMF(E', E) 



poo 

- / Sp(E)RMF(E , ,E)dE 
Jo 



(1) 



For Chandra the RMF matrix is not symmetric. RMFs 
for the aimpoints of the main FI and BI chips are shown 



in Fig. [3] The incident photon energies are binned from 
0.1 keV to 11 keV in 10 eV steps. Thus for FI chips 
there are 1090 raw energy channels; for BI chips there 
are only 1078 raw channels because these chips are not 
calibrated below 0.23 keV. The PHA/PI channel number 
is chosen to be 1024 for both kinds of chips. The values 
in a row, corresponding to a single raw channel or real 
photon energy, correspond to the probability of a photon 
with a given energy being detected in a given PHA/PI 
channel. 

The energy to channel conversion is ideally a one to one 
correlation. But the CCD has a finite energy resolution 
(AE ~ 40 — 170 eV, depending on energy, location and 
CCD type). This is larger than the PHA/PI channel 
width (~15 eV). Moreover, due to instrumental effects 
(escape peak, fluorescence peak low energy peak and tail, 
etc.) a photon of a given energy has a finite probability 
of being detected in channels corresponding to a lower 
energy. This is comparable to "red leak" in optical filters. 

Taking a slice of the RMF for a given energy (7.5 keV 
in Fig. [5]) results in an approximately Gaussian shaped 
curve for the main photo peak, containing ^95% of the 
total, with the FWHM of the main peak contributing 
about 76% of the total. At around channel 20 there 
is the so called low energy peak, with an amplitude of 
about 3 orders of magnitude less than the main peak, 
that contains about 0.1% of the total. If the photon 
energy is high enough, i.e. above the Si K edge, there 
is a secondary peak, the Si fluorescence peak, centered 
on channel 119 (1.73 keV). The maximum of this sec- 
ondary peak is about 2-3 orders of magnitude smaller 
than the main peak. This peak also contains about 0.1% 
of the total. For photon energies above ~2 keV there 
is a third peak, the Si escape peak, with an amplitude 
two orders of magnitude smaller than the main peak. It 
contributes about 0.6%. It follows about 100 channels 
behind the main photo peak. Starting at lower energies 
from the main peak and between the other peaks there 
is a low energy tail at a level about 3-4 orders of mag- 
nitude below the main peak. This tail contains ~4.3% 
with the biggest contribution, ~3%, in between the Si es- 
cape peak and the photopeak. The numbers are almost 
identical for the two kinds of chips. All these features, 
peaks and tail, are due to detector effects. In particular, 
the low energy tail and peak depend on the location at 
which an X-ray photon interacts with the material of the 
CCD and produces an electron cloud. An electron cloud 
from an X-ray photon interaction completely within the 
gate insulator of the CCD chip produces the low energy 
peak. Electron clouds that extend partly into the actual 
detector material, the depleted silicon, contribute to the 
low energy tail. Finally, the main photo peak is made 
up of X-ray photon interactions completely within the 
depleted silicon. The Si fluorescence peak is produced 
by X-rays exciting Si K shell electrons in the detector 
material. The escape peak is produced by fluorescence 
photons that leave the Si substrate or interact away at 
another location with the detector. Fig. 0] shows a slice 
of the RMF for a photon of 7.5 keV for the FI ACIS- 
13 chip (solid line) and for the BI ACIS-S3 chip (dotted 
line) . 

For calibration purposes all Chandra ACIS chips are 
divided into tiles of rectangular form and three different 
sizes: ACIS-S3 (BI) has a tile size of of 32x32 pixels, 
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Fig. 3. — Response matrices for ACIS-13 (top) and ACIS-S3 
(bottom) at the aimpoint. Grey-scale is logarithmic. Numbers 
on the right are energy in units of keV, numbers at the top are 
PHA/PI channel numbers. The vertical dash-dotted lines delimit 
the energy bands. 

ACIS-S1 (BI) has 64x64 pixels, and all the FI chips have 
32x256 pixels. This amounts to 2304 tiles over all chips, 
and thus, in principle 2304 possibly different RMFs. A 
schematic of the tiling is shown in Fig. [T]in the Appendix. 

The main variation among the ACIS RMFs is the 
widening of the main Gaussian peaks with increasing 
distance from the readout for the FI chips, i.e. the en- 
ergy resolution decreases the further a source is from the 
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Fig. 4. — Slice of the response matrices in PHA space for a photon 
of energy 7.5 keV. ACIS-13 (FI) is the solid line and ACIS-S3 (BI) is 
the dotted line. Both curves are for the aimpoint. The photopeak 
is the primary signal; the other features are reduced by some two 
orders of magnitude or more (see text). The broader photopeak 
for ACIS-13 results from the larger distance of the aimpoint to the 
readout compared to ACIS-S3. 

readout due to charge transfer inefficiency. The energy 
resolution changes from a FWHM of ~4 channels (^60 
eV) at 1 keV at the readout to ^9 channels (~130 eV) at 
the opposite side of the chips. For BI chips spatial vari- 
ation of the energy resolution has a more complicated 
shape, but the variation is not as strong as for FI chips. 
At 1 keV the FWHM changes from 6.3 channels at the 
readout to 7.4 channels at the opposite side of the chip. 

For more details a bout the Chandr a RMFs see 
iPrigozhin etaD (|1998f ): iBautz et all (| 1999ft . 

2.4. Photometric bands 

The equivalent of filter shapes for the energy bands in 
Table [T] can be obtained by the convolution of the effec- 
tive area (ARF) with the RMF PI channels over each 
energy band. The filter bands constructed for the aim- 
points of ACIS-13 and ACIS-S3 for 1999 and 2007 are 
shown in Fig. [5] The units are in cm 2 . The four bands 
have quite different normalizations. The soft band, like 
the U band in the optical, has the smallest area. The 
medium-hard and hard bands are quite symmetric and 
have sharp edges to their main response, as is desirable. 
The soft and first medium filters have peaks biased to- 
wards the high energy ends of their bands. This is the 
result of the strong energy dependence of the effective 
area. The low energy leak is quite small. The leakage 
between energy bands is only in the few per cent range, 
and moreover can be corrected quite easily, see below. 

The X-ray filters have a well defined spatial (from 
RMF) and temporal (from ARF) dependence. The spa- 
tial dependence is related to energy resolution, whereas 
the main temporal effect is decreasing sensitivity due to 
accumulation of material on the CCDs. 

Except for the hard band, there is a strong effect due to 
the increased absorption on the CCD between year 1999 
and 2007. The increasing energy resolution towards the 
readout leads to a decreasing overlap between the filter 
bands. However, the energy resolution of Chandra is gen- 
erally sufficiently good that this is a small effect (~3%) 
even at the opposite chip end from the readout. 

These Chandra X-ray filters differ from optical filters in 
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Fig. 5. — X-ray filters for the passbands defined in Sec. I2.ll 
(dotted vertical lines). Filters are constructed for the aimpoints of 
ACIS-13 for 1999 and 2007 (top) and ACIS-S3 for the same years 
(bottom). Dashed curves refer to 1999, solid curves to 2007. 

that they have an area normalization given by the ARF, 
but they are comparable in basic shape. Fig. [5] shows 
a comparison between the X-ray Chandra filters at the 
aimpoint of ACIS-S3 in 2008 (in keV) and SDSS ugri 
filters for zero airmass (in eV). Based on the fractional 
width of the filters, the X-ray filters are only a factor of 
2-4 wider than the SDSS filters. Their variation of the 
peak sensitivity is about the same. 

Due to the effects described in Sec. 12. 31 and illustrated 
in Fig. [4]each filter contributes somewhat to its adjacent 
filters and to all lower energy filters. Thus the filter areas 
are described with a n x n matrix where n is the num- 
ber of filters. Since no photon can be detected at higher 
energies than its own (plus energy resolution effects) the 
upper right elements of this matrix are zero. This ig- 
nores pileup effects in which two photons are detected 
as one at the sum of the individual photon energies. As 
an illustration we give the area matrix at the ACIS-S3 
aimpoint for 2008: 

/18.7 2.1 0.0 0.0 \ 
. _ 4.8 241.1 11.9 0.0 , s 

A ~ 1.2 13.4 998.0 9.9 [ > 

\ 2.1 6.3 30.2 2376.2/ 

The units for the elements are in keV cm 2 . 

3. SPECTRAL MODELS 

X-ray sources can be separated in various categories, 
starting with the distinction between point sources and 
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Fig. 6. — Bottom: Filters for the Chandra X-ray photometry 
system, i.e. combination of ARF and RMF. Abscissa in keV (log 
scale). Top: SDSS ugri filters for airmass. Abscissa in eV (log 
scale) . The variation of peak transmission among the X-ray filters 
is comparable to the variation among the optical filters, with the 
exception of the softest X-ray band. Measured as AE/E the X- 
ray bands are only a factor of few broader than the optical bands. 
Note that the width of the four optical filters covers about 3 eV, 
whereas the four X-ray filters cover about 7000 eV. 

extended sources. The definition of point source is obvi- 
ously a question of angular resolution, e.g. we know that 
in X-ray binaries, which are point-like even with Chan- 
dra , there are sometimes various regions contributing 
to X-ray emission. However, given the excellent angu- 
lar resolution of Chandra (0.3" FWHM at the aimpoint), 
we consider that a point source for Chandra will be a 
point source in 15-20 years as well. The approach of 
photometry is of limited usefulness for extended sources, 
i.e. plasmas, since a definition of regions is somewhat ar- 
bitrary and interpreting photometric differences among 
regions of one source would require plasma diagnostics 
that are beyond the scope of this work. Since plasmas 
are also present in "point sources", e.g. extragalactic 
SNRs, a model is nevertheless included in the following. 

To make the photometry system useful for as large a 
part of the community as possible, it is necessary to study 
and characterize the accuracy of the photometry for a 
wide variety of commonly encountered spectral models. 
Moreover, one of the motivations for a photometric sys- 
tem in X-rays is to not make assumptions about spec- 
tral models. However, the conversion of counts, effec- 
tive area, and exposure time to a flux estimate requires 
some knowledge of the behavior of spectra. We thus 
start by taking common spectral models and define a 
wide range of interest for their parameters. In effect 
this wide ranging combination of spectral shapes and 
parameters represents our ignorance of the true spectral 
shape/parameters of a given source. 

The spectral parameters and models are quite simple. 
For example a flat or highly absorbed power law will 
produce most detected photons in the hard band, while a 
steep or weakly absorbed power law will give a spectrum 
dominated by the soft band. Spectral shapes outside 
the ranges used are extremely hard to detect with X-ray 
telescopes operating in the 0.1-10 keV range. But as we 
will show below the resulting correction to the count to 
flux conversion is quite robust against differences in the 
spectral model in a given energy band. 
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TABLE 2 

Spectral models and parameter ranges used to estimate 
photometric accuracy 



Component 


Parameter 


Range 


Step size 


Absorption 






0.2 in log 


Power law 


r 


-1.0-4.0 


0.2 


Black Body 


kT* 


0.1-2.0 


0.2 


Bremsstrahlung 


kT 6 


0.5-6.0 


0.3 


Opt. thin plasma 


kT 6 


0.1-5.0 


0.3 




abundance 


0.1-1.0 


0.1 


Gaussian Line 


Energies 6 


0.8, 1.3, 1.85, 6.4 


n/a 




Line widths 6 


0.1, 0.5, 1.0 


n/a 



a — units in cm -2 

b - units in keV 

b - units in solar abundance 

4. FLUX ESTIMATION 

In an ideal case the flux F of a source could be sim- 
ply computed from the number of observed counts C s , 
effective area A e f f , and exposure time T as 

F = /V (3) 

A eff T 

However this is true only for flat effective area and no 
cross-talk between energy bands. To obtain a more ac- 
curate estimate of the flux we have to correct for these 
effects. The effect of cross-talk between bands can be 
estimated to first order by computing the contributions 
of individual filters to other energy bands; in reality con- 
tributions to adjacent bands are the dominating factor, 
see Eq. [5] The effect of varying effective area in a band 
is coupled with the effect of the unknown spectral shape 
and is estimated as follows. 

4.1. Correction factors 

Having a spectrum, ARF, and RMF is sufficient to 
compute a theoretical correction factor. This is not the 
general conversion factor from counts to flux that we ulti- 
mately want, but is rather a theoretical value describing 
how accurately finite energy bands can describe a true 
flux for a given spectrum. The flux density at an ob- 
served energy E 1 for a known spectral shape, ARF and 
RMF is 

T(E') = Sp(E) ARF(E) $ RMF(E', E) 

= / Sp{E) ARF(E) RMF(E', E)dE, (4) 
Jo 

where Sp is the original incident spectrum in units of 
photons/cm 2 /s/keV, ARF the effective area in units of 
cm 2 , RMF is the energy-to-channel conversion, which is 
dimensionless. 

Using a broad photometric band we can estimate the 
number of counts in that band C s as the product of the 
integrals over the spectrum and the effective area over 
the energy band. 

C s = f 2 Sp{E)dE- ( 2 ARF(E) dE/AE, (5) 

with Ei and E 2 the lower and upper bound of the energy 
band, and AE the width of the energy band. Note that 
without loss of generality exposure time has been set to 
unity. 

A minor complication is that the observed number of 
counts in a band is not exactly the number of incident 



photons in that energy range due to the redistribution 
effect of the RMF. This can be most easily seen in Fig. 
[5] as the overlap at the filter boundaries. Therefore the 
observed number of counts in any band has contributions 
from other bands as well, particularly for the soft bands. 
In practice the strongest effect (~few to ten per cent) 
is due to adjacent bands. To more distant bands the 
contribution is much smaller (~1%). Ideally the matrix 
of contributions should be diagonalized. But some of the 
matrix elements are zero since softer bands contribute 
only to the next higher band due to energy resolution 
effects. Because of this limitation and since only adjacent 
bands are important contributors of counts, we correct 
for this effect in an iterative way starting at the highest 
energy band, assuming that it does not significantly lose 
counts to even higher energies. Starting at the highest 
band we compute the corrected number of counts as: 

Ain 



r f 


-C S j 


-c si 


Cs\ 


= C s i 


+ C s i 




= Cst 




c sj 


= c s ] 


+ C S ] 



A/, 
Aij 
A, i 

A, 
A 



3 » 



(6) 
(7) 
(8) 
(9) 



where C si j are the observed counts in bands i and j with 
j = i — 1, and A^ the ij-element of the filter area matrix 
(see Eq. [2]). If any of the corrected counts becomes neg- 
ative the corrected counts are set to zero. This iteration 
can trivially be extended to more than adjacent bands. 

The correction factor for a source with known spectral 
shape in a given energy band i is thus defined as 



Ki = 



Cs 



(10) 



The resulting correction factor is dimensionless. Note 
that this formula does not use the information given by 
adjacent bands, which can be used to make a 'color cor- 
rection', as in optical photometry. We will address this 
possibility in a subsequent paper. 

Finite CCD spectral resolution, represented by the 
RMF, and the separate averaging over spectrum and 
ARF generate deviations from the ideal case. This sepa- 
rate averaging over spectrum and ARF gives values close 
to unity only if both are not strongly variable within the 
band. Thus the correction factor is also a measure of 
how flat spectrum and ARF are within the band. In the 
hard and medium-hard bands this is mostly a reason- 
able approximation, but in the soft, and to some degree 
in the medium-soft band, the deviation from flatness is 
large, see Fig. [2j Especially on the FI chips the ARF 
drops precipitously below about ~0.6 keV, by a factor of 
^30 between 0.5 and 0.3 keV. The corresponding drop 
for the BI chips is only a factor ~2.5. Photoelectric ab- 
sorption in the incident spectrum is also a strong factor 
in introducing steep slopes in the observed spectrum at 
low energies, e.g. absorption of ~ 1.5 • 10 22 cm" 2 results 
in a drop of 17 orders of magnitude in flux from 0.5 keV 
to 0.3 keV. 

The correction factor varies as a function of spectral 
parameters in a well behaved way. As an example Fig. 
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[7| shows the correction factor surface for a power law 
incident spectrum in the four photometric bands at the 
aimpoint of the BI ACIS-S3 chip in 2007 as a function of 
absorption and photon index. 

In general the correction is relatively small (~20-30%). 
The main variation occurs for large absorption values 
(Nh > 10 23 cm -2 ) where the correction factor drops 
to zero in the soft and medium-soft bands, or strongly 
increases in the medium-hard band. For practical ap- 
plication, however, a correction factor of zero is not a 
problem because a zero correction means that there is 
no flux in that energy band down to computational ac- 
curacy (J SpdE — in Eq. fTDj) . thus a correction is not 
useful. Practically, in the code we set incident fluxes to 
zero if the integral over the spectrum in an energy band 
is smaller than 10~ 30 photons/cm 2 /s. 

Importantly, the shapes of the correction factor sur- 
faces in these plots are similar for different spectral 
shapes, see Figs. [21 [31 H]in the Appendix. The difference 
between two correction surfaces in one band is generally 
less than ~15%, and never worse than ^40% for the most 
extreme spectral shapes. Since the spectral shape of a 
source is a priori unknown, and the point of this work is 
to not assume any shape, we combine all the points on 
each grid for all the different spectral shapes. The dis- 
tribution of correction factors taken from all spectra is 
shown in Fig[5]for the aimpoint of the BI ACIS-S3 chip 
in 1999 and 2008. The histograms are strongly peaked, 
although they have significant tails towards larger correc- 
tion factors. These tails are especially pronounced in the 
medium-hard band and are due to spectra with strong 
absorption values. 

Ignoring spectra with large absorption values (Nh > 
10 23 cm~ 2 ) in the computation of the correction factor 
has no significant impact on the peak of these distri- 
butions but considerably reduces the asymmetry of the 
distribution for the medium-hard and hard band; the 
tails almost disappear. This property of highly absorbed 
spectra may be used in a "color correction" similar to 
optical photometry. E.g. if there are no counts in the 
soft and/or medium-soft band, this is already an indica- 
tion of high absorption, which may be used to change the 
correction factor distribution in the medium-soft and/or 
medium-hard and hard band. We address this band ratio 
correction in a subsequent paper. 

Given our ignorance about the intrinsic spectral shape 
of a source the correction factor distribution can be con- 
sidered the probability density function PDF of picking 
the right correction factor for conversion from counts to 
flux. Thus the probability density function for the cor- 
rection factors is simply the conjunction of all correction 
factors Ki for a given spectral shape i: 

PDF(K) = f\ Ki. (11) 

i 

with N sp the number of different spectra, in this case the 
sum of parameter combinations from Table [H and Ki the 
correction factor for spectral parameter combination i 
according to Eq. [TU] Note that N sp has to be sufficiently 
large to obtain a well behaved distribution. With our 
choice of spectral shapes and parameters N sp is 1454. 

It is important to note here that our selection of spec- 
tral shapes serves effectively as a prior on this distribu- 



tion. This prior is flat in the spectral model/parameter 
space, i.e. all model/parameter combinations are given 
equal weight in computing the correction factor distribu- 
tion. If the correct spectral shape is not covered by this 
range the method is no longer valid, which is the reason 
for the large parameter space we use. 

Considering the correction factor distribution as a PDF 

we use the mode of the distribution PDF as the cor- 
rection factor. As an error estimate we use the range 
that covers 68% of the distribution as measured from 
the mode. The mode of the correction factor distribu- 
tion with the 68% confidence level (CL) errors at the 
aimpoint of ACIS-S3 for each band for the past 10 years 
is shown in Fig. [9] 

If not noted otherwise all error bars and quoted errors 
in the following correspond to the 68% confidence level. 

4.2. Flux 

Using the knowledge of the effects of separate averag- 
ing over spectrum and ARE, our ignorance of the true 
spectral shape of a source, and our estimates for correct- 
ing these effects, we can compute a flux estimate F for a 
source in energy band i at a given location on ACIS as 

Ft = ^ (12) 

/ (ARF <g> RMF)idE ■ T ■ PDFi 

with C' s as the background subtracted counts corrected 
for filter overlap in energy band i, ARF and RMF the ef- 
fective area and redistribution matrix in band i, T the ex- 
posure time, and PDFi the mode of the correction factor 
distribution in band i. Fi is in units of photons/s/cm 2 . 

4.3. Time variation 

Theoretically for Chandra all points on the chips have 
their own value for the correction factor which in addi- 
tion is time dependent. Practically, however, the spatial 
variation across a chip is much smaller than other effects. 
E.g. the lack of knowledge of the correct spectral shape 
which is encoded in the variance of the correction factor 
distribution is of order 15%-20%. Whereas the variation 
due to changing energy resolution across a chip is about 
1%. 

Instead temporal changes are important. The change 
in sensitivity with time significantly affects the effective 
area and thus the correction factor. Fig. [9] shows the 
change in the correction factor at the aimpoint of ACIS- 
S3 for the different energy bands between 1999 and 2007. 
There is significant change in the correction factor for the 
soft band up to about 2002 due to accumulation of ma- 
terial. The total change is about 28% between 1999 and 
2008, however the medium- hard and hard band are virtu- 
ally unaffected. The variations for the ACIS-I3 aimpoint 
are almost identical, except the variation in the soft band 
is only about 20%. 

5. UNCERTAINTIES 

The uncertainties in the X-ray photometry can be sep- 
arated into statistical and systematic errors. The sta- 
tistical errors are due only to counting statistics, and so 
are observation dependent. Due to the fact that our pho- 
tometry method also deals with small number of counts, 
Gaussian statistics and error propagation are not appro- 
priate tools. Independent of counting statistics, we con- 
sider the distribution of correction factors for a specific 
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FlG. 7. — Surface of correction factors versus spectral parameters 
FI ACIS-S3 chip in 2008. 
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Fig. 8. — Distribution of correction factors taken from all spectra 
for aimpoint of ACIS-S3 in 1999 and 2008. Soft band-solid his- 
togram, medium-soft band-dash-dotted histogram, medium-hard 
band-dashed histogram, hard band— dotted histogram. 

chip location and time as a systematic error. This distri- 
bution is dependent on the number and kind of spectral 
shapes and the chip location in time and space. 

5.1. Systematic uncertainties 

Neither for ARF nor RMF values uncertainties are pro- 
vided although they are estimated to be below 10%, and 
around 5% for mo st of the energy range for the ARF 
(|Drake et al.ll2006h . We therefore use the correction fac- 
tor distribution to estimate systematic uncertainties in 
our program. Individual correction factors are the result 
of combining spectral shape, ARF, and RMF. Since we 
use a wide range of spectral shapes and parameters the 
influence of the spectral shapes/parameters on the uncer- 
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Fig. 9. — Change of correction factor with 68% error bars at 
the aimpoint of ACIS-S3 between 1999 and 2008. The top panel 
is the hard band correction factor, second from top medium-hard 
band correction, third from top medium-soft band correction, and 
bottom panel the soft band correction. For the soft band the mode 
of the correction factor distribution is at the lowest value, thus 
there are no lower error bars (see Fig. [§}. 

tainty in flux is certainly much larger than uncertainties 
in the ARF or RMF. 

To obtain the systematic error we choose the mode of 
the correction factor distribution at a given location of 
the instrument as the correction factor for that location 
and integrate from that position in positive and nega- 
tive direction until 68% of the distribution are covered. 
In general this will result in asymmetric errors since the 
distribution is skewed. The flux is inversely proportional 
to the correction factor, thus the systematic error does 
not have to be propagated beyond inverting and multi- 
plying by a numerical factor. 

5.2. Statistical uncertainties 

To obtain the uncertainty on the flux in a band we 
propagate the errors on the source counts through the 
corrections applied to the observed counts in a band. In 
th e presence of backg round we use the method proposed 
bv lKraft et alJ (|1991l ) to compute the number of source 
counts and the uncertainties in that number. Note that 
the uncertainties are asymmetric and not exactly Pois- 
sonian. 

Unfortunately there is no obvious or generally accepted 
way to prop agate asymme tric errors. Here we follow the 
approach of lBarlowl (|2004f ) for combining asymmetric er- 
rors. The approach is based on the idea to parametrize 
the Log-likelihood curve of the original probability dis- 
tribution with 3 parameters: Location, scale, and skew 
for each measurement. Then the Log- likelihood func- 
tions for individual measurements are combined, in this 
case added. Obviously there are numerous ways to 
parametrize a function through three points. In our pro- 
gram we use the parametrization described as Gaussian 
with linear variance. This choice is purely empirical and 
made because this parametrization approximates very 
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well a Poissonian distribution. And although the back- 
ground subtracted counts and errors are not Poissonian, 
the difference is quite small. The likelihood function for 
each part is then given by 

/(a) = ln(L(x; a)) = -~ -f (13) 

2 <J-(7+ + (<7+ — cr_ja 

where a is the value, and cr_) are the asymmetric uncer- 
tainties in negative and positive direction. 

Combining the likelihood functions for two distribu- 
tions one can then obtain the combined uncertainties at 
the locations where / = /(d) — |, with a being the sum 
of the counts. Using this method we obtain uncertain- 
ties that are slightly larger than assu ming Gaussian error 
propagation for the Gehrels errors (Gehrels 1986). For 
large number of counts the error approximates the Gaus- 
sian expectation value. Thus we consider our statistical 
errors to be conservative. 

6. VALIDATION 

We compare the results of our program with simulated 
and real data to validate the results and check for varia- 
tions. 

6.1. Simulated data 

We simulate 100 spectra each for parameter combina- 
tions of two spectral shapes, a black body with a tem- 
perature of 0.3 keV and and power law with a photon 
index of 2. Both spectral shapes have absorption column 
densities of 10 20 and 10 22 cm" 2 . All spectra contain 1000 
counts and no background. This is a good approximation 
for Chandra point sources. The spectra are fitted with 
Xspec and fluxes are computed using the flux command. 
These are compared with results from our method. 

Results from this comparison are shown in Fig. 1101 
The figures show Xspec fluxes (F/it) over photometric 
fluxes (-Fp/jot) computed with our method versus photo- 
metric fluxes for the two spectral shapes and four spec- 
tral parameter combinations. Different bands are labeled 
and shown with different symbols. The error bars used 
correspond to 68% CL. Ideally all y-axis values would be 
one. 

It is clear from the figures that this photometry com- 
putes fluxes within better than 30% in most cases, even 
in the absence of band ratio (color) corrections. That the 
medium-hard band (M2) has the best correspondence be- 
tween X-ray photometric fluxes and Xspec fluxes is sim- 
ply based on the fact that the medium-hard filter has 
the best properties, very steep edges and a relatively flat 
top (variation of ^10%) as shown in Fig. [5l Even in the 
worst cases it is better than 50%. It should be noted that 
even in the worst cases the systematic error of 68% CL 
covers the discrepancy between the true and estimated 
fluxes. The large discrepancies are due to several rea- 
sons: (1) a relatively strongly varying spectrum in the 
energy band, (2) a varying ARF, and (3) differences be- 
tween the correction factor for that spectrum and the 
mode of the correction factor distribution. The last fac- 
tor is quite small, at least for the spectral shapes chosen 
above. The main contribution to the discrepancy comes 
from the combination of varying spectrum and ARF and 
a high number of counts. 

With a high number of counts (>^200) in a band, 
the changes in the ARF over an energy band become 



important, particularly for the broad hard band, and the 
resulting flux differences become significant with respect 
to the statistical error. 

6.2. Real data 

The real data we use for our comparison are relatively 
bright (> 100 source counts) sources from Chandra ob- 
servations of M 33. We use two approaches to estimate 
the flux independently from our method. First we fit 
the source spectra with spectral models and use the flux 
command in Xspec to g e t the photon flux; for best fit 
models see iGrimm "eTall (|200l . 

Secondly we take each photon in an energy band, use 
the inverse of the ARF at the photon energy, sum this in- 
verse over all photons in the energy band, and divide the 
result by the exposure time which gives the photon flux. 
This is the simplest possible photometry method. How- 
ever, it requires knowledge of individual photon energies, 
whereas our method works with integrated quantities. 
This means that the accuracy of this method depends 
not only on the number of counts, but for weak sources 
also on where in the energy band the few photons fall. 
And moreover, it does not take into account effects of 
the RMF which are important for spectra that produce 
very different counts in different bands. 

Results for sources with well fitted spectra are shown 
in Fig. QT] The figure shows the flux estimates of our 
methods as red points, Xspec results as green squares, 
and summation of inverse effective area as blue triangles. 
The error bars are 68% CL statistical errors only. In gen- 
eral there is good agreement between all three methods 
of flux estimation. 

It is apparent from the figure that there are sometimes 
significant discrepancies between the flux estimate from 
our photometry and Xspec fluxes. However the summa- 
tion of the inverse effective area and our method agree 
quite well in general. So it is the spectral fitting that re- 
sults in relatively large deviations from the flux estimates 
obtained by the other methods. This is probably due to 
several factors. First, the spectrum is fitted in Xspec us- 
ing the whole energy range. This can result in an undcr- 
or overestimate of fluxes in smaller energy bands, that do 
not strongly contribute to the overall fit. Furthermore, 
it can be the case that the chosen spectrum is not a good 
representation of the data but due to insufficient quality 
of the data this is not apparent from bad fits. E.g. the 
spectrum used for the left panel in Fig. [IT] has a reduced 
X 2 of 0.4 which clearly indicates that the bremsstrahlung 
spectrum chosen for the fit is not the true representation 
of the data. 

7. CONCLUSION 

We have presented a system of X-ray photometry for 
the Chandra satellite. The system we have developed 
relies on a knowledge of effective area and the energy- 
to-channel conversion to construct X-ray filters, but is 
unbiased by assumptions about the spectral shape of a 
source. We have shown that the filters are comparable 
to filters in the optical and infrared, and that our pho- 
tometric system in X-rays is able to estimate fluxes to 
within about 20%. Even in the worst cases it is better 
than 50%. We have incorporated methods to estimate 
systematic errors and consistently propagate statistical 
as well as systematic errors. 
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Fig. 11. — Flux estimates of our X-ray photometry as (red) 
circles, Xspec results are (green) squares, and summation of inverse 
effective area as (blue) triangles versus energy band. Error bars are 
68% CL statistical errors. 

Due to the construction method employed our fil- 
ter system is very flexible and can be adapted readily 
to other CCD X-ray detectors, in particular to XMM- 
Newton EPIC. The code to compute fluxes is available 
at a web page 2 . A table with a selection of correction 
factors for ACIS-S3 is given in Appendix A in the elec- 
tronic edition and available at the same URL. The table 
contains only every fourth tile because of the generally 



slowly varying correction factor with chip location. In 
the future we will explore potential improvements to X- 
ray photometry by means of: 

• Making color corrections using band ratios. A pre- 
liminary investigation suggests that extreme spec- 
tra (e.g. highly absorbed or high photon indexes) 
will gain significantly in the accuracy of flux esti- 
mates and even normal spectra will have a reduced 
error range. 

• Optimizing the choice of bands. The properties of 
our current method show that there is a correlation 
between the accuracy of flux estimates and the fil- 
ter shape. The more "boxy" a filter is the better 
the flux estimate will be. This suggests a limit for 
the width of a filter at which point deviations from 
boxiness result in an accuracy of the flux estimate 
below a certain value. 



• Comparing results for Chandra ACIS and XMM- 
Newton EPIC data. Chandra and XMM have simi- 
lar instrumental setups (ARF and RMF) and over- 
lapping science capabilities. And given the large 
numbers of sources observed by these two X-ray 
missions it is very important to be able to compare 
results for the two. 
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APPENDIX 




Fig. 1. — Schematic of the tiling of the full ACIS detector. Direction of the x-axis is indicated by the arrows. The y-axis is counterclockwise 
to the x-axis. Aimpoints for the 13 and S3 chips are shown as crosses. The readout for chips is on the side of the arrows. 



X-ray photometry for Chandra 
TABLE 1 Correction factors for selected tiles on ACIS-S3. 
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a - X and Y arc the tile center in chip coordinates. 

b - Values arc the mode of the correction factor distribution (EgfTT) with the 68% CL error (sec Sec. f4"Tl fc FTP ■ The values should be used as input in Eq. ^] 
c - Aimpoint on ACIS-S3 
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